A quantum computer using a trapped-ion spin molecule and microwave radiation 
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We propose a new design for a quantum information processor where qubits are encoded into 
Hyperfine states of ions held in a linear array of individually tailored linear microtraps and 
sitting in a spatially varying magnetic field. The magnetic field gradient introduces spatially 
dependent qubit transition frequencies and a type of spin-spin interaction between qubits. Single 
and multi-qubit manipulation is achieved via resonant microwave pulses as in liquid-NMR quantum 
computation while the qubit readout and reset is achieved through trapped-ion fluorescence 
shelving techniques. By adjusting the microtrap configurations we can tailor, in hardware, the 
qubit resonance frequencies and coupling strengths. We show the system possesses a side-band 
transition structure which does not scale with the size of the processor allowing scalable frequency 
discrimination between qubits. By using large magnetic field gradients, one can reset individual 
qubits in the ion chain via frequency selective optical pulses to implement quantum error correction 
thus avoiding the need for many tightly focused laser beams. 
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A quantum computer requires well-characterized 
quantum bits (qubits), sufficient control to perform 
a universal set of quantum gates, long decoherence 
times compared to gate times and a means of per- 
forming qubit readout To be of practical use, 
any design must scale to hundreds of qubits without 
enormous technical overhead. Current ion trap and 
nuclear magnetic resonance (NMR), implementations 
are capable of satisfying most of these criteria to varying 
degrees. Liquid-state NMR quantum computing has 
demonstrated very precise control in the manipulation 
of qubits with microwave and radio frequency pulses and 
the execution of a number of quantum algorithms. The 
main drawback of Liquid-NMR quantum computing is 
the extreme difficulty, due to tiny spin polarization and 
weak measurements, in scaling up to a large number of 
qubits. Ion trap quantum computing has recently shown 
progress in implementing a quantum algorithm 0, and 
a controlled interaction between two qubits Q. Here 
the initialization and readout of qubits is particularly 
good. However, there are barriers to be overcome such 
as decoherence due to ion heating, very precise laser 
focussing and stability problems of operating at optical 
frequencies. These difficulties may hinder the ability 
to scale optically-addressed ion-trap technologies to 
hundreds of ions although a number of potentially scal- 
able designs involving moving and stationary optically 
manipulated trapped ions have been proposed Q. In 
this article we describe a new design for a potentially 
highly scalable quantum information processor which 
combines the trapped-ion and NMR technologies in a 
manner that retains the advantages of both. 
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We assume that two Hyperfine levels in each ion serve 
as the qubit. As demonstrated by Wunderlich et. al. 
HH0, one can induce a coupling between the qubits 
through the application of a magnetic field gradient 
along a string of ions. This, together with the Coloumb 
force, couples the qubits together in a way analogous to 
the spin-spin coupling observed between nuclear spins 
in a molecule. Obviously then, similar techniques to 
those used in NMR, and the related field of electron 
spin resonance (ESR), can be utilised to realize two 
and multi-qubit quantum gates. The work presented by 
Wunderlich et. at, though possessing obvious benefits 
for the precise manipulation of quantum information, 
however displays a number of serious obstacles with 
respect to scalability. In that design we will find that the 
spin-spin coupling varies in strength throughout the ion 
chain thus leading to varying quantum gate durations 
between qubits depending on their location in the ion 
chain. This, itself poses a serious difficulty in that a 
scalable device will require device- wide parallel quantum 
information processing and quantum error correction 
and this will need careful synchronisation of logical 
operations throughout the device. More problematic 
is the existence in the Wunderlich design of an upper 
limit to the number of qubits which can be cleanly 
frequency differentiated for a given size of magnetic field 
gradient. Finally, the technical difficulty and potential 
added source of decoherence introduced by the focusing 
of individual laser beams on each ion in the chain for 
individual qubit initialization, readout and reset may 
also be problematic for scalability 0. Here we show 
that by considering a design where the ions are stored 
in a linear array of individually tailored linear micro- 
traps 0, we can surmount all of the above mentioned 
difficulties and achieve significant gains in scalability. 



2 



In particular we find that the resulting system is, as 
before, analogous to an N atom molecule with spin-spin 
couplings, but now these couplings can be tailored 
very precisely, and quite robustly, through altering the 
individual microtrap parameters. The new design is also 
exactly analogous to a true ion crystal where each ion 
is harmonically bound to a periodic spatial lattice. The 
resulting qubit resonance frequencies together with the 
vibrational side-transitions form a band structure (as 
in condensed matter systems), where a band's extent 
now only depends on the inter-ion spacing and end trap 
strength and does not grow with the number of ions in 
the chain. Finally, as previously shown by Wunderlich 
@, the magnetic field gradient allows the qubits to 
be individually frequency-addressed in the microwave 
via Zeeman splitting of the Hyperfine structure. We 
further show that for large magnetic field gradients one 
can engineer for frequency-addressed optical resetting 
of an individual Hyperfine qubit via Zeeman splitting 
of the optical transitions. This capability allows for 
scalable quantum error correction and relaxes the need 
for focusing of individual laser beams on each ion 
through the use of frequency multiplexed laser sources. 
Previous difficulties predicted in liquid NMR quantum 
information processing related to the saturation of the 
available RF-bandwidth with a moderate number of 
frequency differentiated qubits is no longer an issue here 
as the qubit manipulation in the MW, (or readout in 
the optical), has available radiation sources which are 
tunable over many MHz (GHz). Further, our design 
uses the highly developed and relatively widespread 
technology of microwave pulse synthesis to precisely 
manipulate the qubits. This avoids the need to develop 
and maintain ultra-stable narrow-band laser systems as 
used in in 0, y( > although high frequency optical systems 
will still be required for the qubit initialization, reset 
and readout processes. Our model differs fundamentally 
from other quantum computer designs based on trapped 
ions. In jlOj , a two-qubit gate is realised through the 
exchange of a phonon through collective vibrational 
motion of the ions while in the ions are also stored 
in an array of microtraps and qubit gates are realised 
through a state-dependent displacement of the motional 
wave-packet of the ions. In the scheme of Wunderlich 
and in our scheme, the spin-spin coupling is achieved 
via virtual excitation of the ion's motion and thus gate 
operations should be more robust against decoherence 
arising from the motion of the ions as in the schemes 
presented in 0, 0, 0] . 

We first explain qualitatively how the model operates 
and follow this with a quantitative analysis. We consider 
N ions, with each ion occupying a separate harmonic 
oscillator potential well, arranged in a linear array. We 
assume we can independently fix (1), the strength of the 
each potential well and, (2) the separation between each 
well. Each ion stores a qubit in two Hyperfine states 
and taken together constitutes a quantum register. The 



initialization of the qubits is performed as usual via 
fluorescence shelving and repumping [12|. A magnetic 
field gradient is applied along the ion string and results 
in frequency differentiated qubits in the string due 
to the Zeeman shifting of each ion's Hyperfine levels. 
Single-qubit operations on a given ion are performed by 
illuminating the ion string with a pulse of microwave 
radiation of the appropriate frequency. The spin-spin 
coupling between the ions is achieved through a combi- 
nation of the Coloumb force and the qubit dependent 
Zeeman energy which is spatially modulated by the 
magnetic field gradient The resulting spin-spin 

coupling is thus "on" all the time and gates are achieved 
via a combination of single-qubit operations and free 
evolution as in liquid-NMR quantum computation. 
Thus one can use the considerable knowledge of NMR 
refocusing and averaging techniques |l6[ . to tailor the 
system Hamiltonian and perform quantum information 
processing. However, quantum computing using liquid 
state NMR and an ion spin-molecule sitting in a single 
trapping potential as proposed by Wunderlich, differ in 
that the latter allows one to control (although in a very 
rough manner), in hardware, the inter-qubit coupling 
strengths. For the latter the resulting inter-qubit 
couplings vary in strength throughout the ion-chain no 
matter how one alters the single trap parameters (see 
Fig. nj . This lack of homogeneity represents a significant 
barrier for the scalability of 6]. In addition, using a 
single linear trap to contain the ion string yields limits 
on the number of ions in the string that can be cleanly 
frequency differentiated for a given size of magnetic 
field gradient. By introducing linear microtraps and 
individually tailoring their strengths and locations we 
can surmount both of the above problems to yield a 
design with high scalability. We can precisely engineer 
the inter-qubit couplings to be homogenous throughout 
the chain while maintaining uniform separation between 
the ions and consequently, well differentiated qubit 
resonant frequencies irrespective of the size of the 
device. Having an homogeneous system enormously 
boosts the scalability of device as the alternative would 
imply that all gate operations (pulse sequences), would 
be highly dependent on the location of the ions in the 
chain which are involved in the gate. In addition by 
varying the ion-separation and microtrap strengths 
we also have some control over the relative sizes of 
non-nearest-neighbour to nearest-neighbour spin-spin 
coupling strengths. In [(|, and in our model, the ions 
in the string are spatially separated on length scales 
varying from 2 — 10/im, and although the individual 
focusing of lasers on each ion is possible, such a readout 
scheme may yield too much decoherence due to laser 
pointing fluctuations, amplitude noise etc. Also in [||, 
one can readout the entire qubit register at once at the 
completion of a quantum algorithm. However to execute 
quantum error correction, intermediate measurements 
are usually desired. However, as shown in 9], quantum 
error correction only requires the capability of supplying 
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fresh ancilla qubits, or in the present case, the resetting 
of any individual qubit to a preset quantum state. 
We show that in the case of large constant magnetic 
field gradients, the resulting spatial dependence of the 
energies of both the qubit Hyperfine levels and the 
optically-excited "readout" level allows one to achieved 
optical frequency differentiated reset through optical 
pumping or through engineered decoherence processes, 
with little disturbance of neighboring qubits. This yields 
another boost for the scalability of the resulting design. 

We now look at the model more quantitatively. Con- 
sider N ions each of mass m confined in N individual 
harmonic potential wells. The i th linear trap has fre- 
quency iOi and is located at position ki along the x-axis. 
A spatially varying magnetic field, B = (B + bx)z, is ap- 
plied across the line of ions. The resulting Hamiltonian, 
to second order in the ion's vibrational motion, is given 

by 0, 



the Hessian A. Firstly the off-diagonal elements of A are 
functions only of the inter-ion spatial separations. We 
can fix the positions of the individual microtraps, ki, to 
achieve a uniform ion separation, h, thus enormously sim- 
plifying the structure of A. The diagonal elements of A 
are functions of both the inter-ion separation, h, and the 
individual trap strengths. Letting = mujf, be a mea- 
sure of the i'th trap strength, ei = (e 2 / Aneoh) / '(^gih 2 ) 
fll| . and assuming only nearest neighbor coupling, one 
can analytically show that for uniform off-diagonal J's 
one must have 
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The first term represents the electronic Hamiltonian 
of the qubits, now with separated qubit resonances 
w n (xo,n), and where xo,n is the equilibrium location 
of the ro'th ion in the string sitting in the trapping 
potential. The second term describes the collective 
quantized vibrational motion of the ions. Even though 
the ions are now in individual linear traps, this term 
describes the quadratic interactions between the ions 
and is essentially identical to the case of a single 
harmonic trap. The last term expresses the pairwise 
coupling between qubits analogous to the well-known 
spin-spin coupling in molecules used for NMR quantum 
computing. In , it is then proposed that this last term 
can be used to implement quantum gates. 

One can show that the spin-spin coupling between ions 
n and m can be expressed in the relatively simple form: 
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where ^ | , is the gradient of the qubit transition reso- 
nant frequency for ion k, at equilibrium location xo,fc; A 
is the Hessian of the potential in which the ions sit; D 
is the unitary transformation matrix that diagonalises A 



and 



are the eigenvalues of A. 



In the high-field or Paschen-Bach limit, (Bq ~ IT), 
the frequency gradients are independent of x and thus n, 
and so J oc A^ 1 . To achieve uniform off-diagonal values 
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with all interior traps being weaker in strength than the 
end traps. Resorting to numerics to include non-nearest 
neighbor terms in the Hessian, one can determine the 
interior microtrap strengths for any given length of ion 
string with a uniform spatial ion separation. As one 
might expect both ki and gi are symmetrically valued 
about the mid-point of the ion chain. As an example, for 
Yb + ions, where the qubit is encoded in the Hyperfine 
states |0) ee |6 Si, M/ = \, Mj = |1) = 

6 Si, M[ — \Mj = i), of each ion we can obtain the 
microtrap strengths in Fig. [21 Once suitable values 
for gi are obtained, the resulting microtrap positions 
ki, which yield the specified uniform ion separation, 
can be separately computed numerically. As expected, 
the microtrap separation decreases slightly towards the 
edges of the ion string to compensate for the larger 
outwards Coloumb forces experienced by ions in these 
regions. From Fig. ^ the resulting improvement in the 
uniformity of the J-couplings throughout the ion string 
is quite striking. In addition however, we have also suc- 
ceeded in significantly reducing the non-nearest neighbor 
J-couplings, thus reducing these potential sources of 
error in the performance of nearest neighbor quantum 
gates. Further, by varying both the inter-ion separation 
h, and end-trap strength g\, we can control the degree 
of suppression of non-nearest neighbor couplings. Thus, 
through tailoring the microtrap configuration we can 
achieve uniform off-diagonal coupling strengths through- 
out the ion-string. Other types of J-couplings may also 
be tailored for. 

Of course it is important to estimate the degree of in- 
tolerance of the design due to misalignments in the trap 
strengths and positions. In numerical simulations, fluc- 
tuations of the ion positions have minimal effect on the 
uniformity of the J-couplings. Given 10 ions in a trap, 
with vx = 1MHz, b = lOOOT/m and h =10/im, setting 



X0.7 



for J ~ A , i.e. J^. 



J 



we now tailor the values of 



Xq n + Sx n , where Xq „ are the microtrap loca- 
tions giving a uniform J-coupling, and 8x n randomly dis- 
tributed deviations within a range ±0.1/im, the nearest- 
neighbour J-couplings vary by approximately 1% on av- 
erage. Considering the spatial extent of the lowest axial 
vibrational mode of each ion is ~2nm, the J-couplings 
maintain their uniformity when heating occurs. For 



4 



variations of the individual linear microtrap strengths, 
v,i = vf + Sfi, with 5vi randomly distributed deviations 
within a range ±10kHz (~ ±1% of vi), results in 1% 
variations of the nearest-neighbour J— couplings. Fur- 
ther, one can employ so-called robust pulse sequences to 
overcome small inhomogeneities in the coupling strengths 
throughout the chain |17| . 

We now take some time to outline the complete 
scalability of the qubit frequency discrimination within 
our design. As noted in [y, [3 , the vibrational energy of 
the highest collective motional mode of N ions in a single 
harmonic trap subject to a magnetic gradient of size 
dB/dx = 6, increases linearly with N, i.e. v max ~ v*N. 
For clean frequency discrimination between two neigh- 
boring qubit resonant frequencies, u>i and Wj+i, one 
requires, — u>i\ oc dB/dx > 2v max ~ 2v*N. 

Thus for a given magnetic field gradient b, there is a 
maximum number of ions, N* , which can be cleanly 
discriminated. In our model, however, the vibrational 
frequencies of the collective motion of N ions, held in 
microtraps tailored for uniform spatial separation and 
qubit coupling strengths J, asymptotes to a given value 
Vmax{N,h,gx) -> v ma x(h,gi), for large N. Empirically, 
defining Av(N,h,gi) = v max (h,gi) - v m ax{N,h,gx), we 
find Av(N,h,gi) ~ Av{Nj,h,g x ) x {N/Nj)' 1 - 78 . 
Empirically, we also can find the scaling of 
this convergence with ion separation h, to be 
Av(N,h, gi )~Av(N,h,gi) x {h/h)-\ 

It is necessary finally, to discuss in more detail how 
it is proposed that the qubits can be individually ma- 
nipulated, read out and reset. First, we will talk about 
manipulation. Since we are proposing to using two Hy- 
perfine levels for the qubit, long wavelength radiation is 
involved and so we will have to manipulate the qubits 
via frequency selection rather than by attempting to fo- 
cus laser light on the ions individually. As an exam- 
ple, we again identify the \Si,Mi = \,Mj = — 5) and 

\Sx,Mi = \,Mj = i), states of 171 Yb+ as the qubit 
|0) and |1). The Zeeman splitting of the Hyperfine lev- 
els lies in the Paschen-Back region when Um3M > 1 

, , &HFS 

where Ehfs is the Hyperfine splitting constant, gj the 
g-factor and B(z) the strength of the magnetic field. For 
Yb + , this gives B > IT. If we are in this region then the 
difference between the manipulation frequencies of neigh- 
bouring qubits is given by Av mampu i ate = gj ^ bSz , where 

b = dB Q^ and 5z is the separation between the ions. 
Choosing a field gradient of 1000 T/m and a typical ion 
separation distance of 10/um we can realise a neighbour- 
ing qubit frequency difference of 280MHz, irrespective 
of the length of the ion-string. This separation is easily 
differentiated by current ESR Microwave pulse spectrom- 
eters. 

The final readout of the state of the qubits can be 
performed, as usual, by spatially resolved resonance fluo- 
rescence. However, scalability will require the capability 
to perform quantum error correction. This capability 



will require individual ion readout and/or reset. How- 
ever, as the spatial separation between ions is ~ fim, 
readout/reset of individual ions through individual laser 
focusing may be difficult. Instead, one can use frequency 
addressing of the optical transitions to achieve individ- 
ual ion readout /reset. For B > IT, both readout lev- 
els, \6Si,Mi = \,Mj = \) and \GPi,Mi = \,Mj = 

— \), are in the Paschen-Back region and the resulting 
frequency separation between the optical transitions of 
neighbouring qubits is Av read out = f ^ Sz ■ This gives 
an optical frequency difference of 187MHz for the above 
situation which is far in excess of the excited state lifetime 
~ 10MHz, and thus cleanly discriminated between neigh- 
boring ions. Quantum error correction typically involves 
the measurement of an individual qubit state but to use 
fluorescence detection in a manner which leaves the spec- 
tator qubits undisturbed requires the additional correc- 
tion of considerable AC Stark shifts suffered by neighbor- 
ing ions during the illumination and acquisition of ~ 10 4 
fluorescent photons from the target ion. Techniques to 
perform this correction have been demonstrated |18| , but 
such invasive detection may not be necessary as only the 
capability to perform a qubit specific reset, are required 
for quantum error correction |ij . Resetting the Hyperfine 
qubit to a known quantum state should only require the 
scattering of far fewer fluorescent photons of photons 
~ M = 10° — 10 1 ), via optical pumping and the result- 
ing phase acquired by a neighboring ion off-resonant by 
L ~ 19, radiative lincwidths, will be <fi ~ M/L ~ 1/5 
radians. The correction of such phases are well within 
the capabilities of the techniques demonstrated in (l8| . 
Alternatively, a technique to reset an individual Hyper- 
fine qubit by simultaneously resonantly irradiating the 
qubit and optical transitions with strong MW and ex- 
tremely weak optical radiation at 369nm (with a laser 
power of nW), has previously been demonstrated in |l9j |. 
This process of engineered decoherence, may also be used 
as a qubit reset. 

Recent applications of narrowband ultraviolet laser 
diodes diode lasers to Ytterbium trapping might be 
adapted for the frequency selective optical addressing 
[20| . By using electro-optic modulators and external 
cavity tuning elements, tuning ranges of several GHz 
are currently possible at these wavelengths. Through 
illumination of the entire ion chain by a frequency 
multiplexed laser source one can achieve individual 
qubit readout and reset in a highly scalable manner. 
Achieving in-trap magnetic field gradients of lOOOT/m 
over several micron can be achieved through methods 
similar to those presented in pl|. Also, engineering 
linear ion traps on the micron scale is possible as 
described in |22j. However, by decreasing h, we increase 
the qubit coupling strengths while decreasing the fre- 
quency separation between the qubit manipulation and 
readout frequencies. As a final example, in Table 1, we 
have computed the parameters for ten Yb + ions, with 
B a = IT, b = lOOOT/m, v x = 1MHz, and h = 10/zm, 
held in individually tailored microtraps for uniform 
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separation and J— couplings. Here the uniform nearest 
neighbour coupling strength, Jj^+i ~ 850Hz. The tab- 
ulated parameters are within current microwave source 
capabilities and clearly indicate scalable frequency 
discrimination for the manipulation, readout and reset 
of neighboring qubits. 

In conclusion, we have presented here a new design 
for a trapped-ion spin-molecule quantum information 
processor which we feel should be highly scalable. Our 
design surmounts the difficulties in scalability faced by 
the scheme of Wunderlich et.al. By storing the ions in 
individually tailored linear microtraps we can engineer 
homogenous coupling strengths throughout the chain 
while simultaneously achieving a constant frequency 
separation between neighboring qubits, irrespective of 



the length of the ion chain. This enormously simplifies 
the microwave pulse sequences required to manipulate 
the quantum information in the ion chain. At the same 
time we have removed the limit on ion numbers which 
can now be cleanly frequency discriminated in a fully 
scalable manner. Finally, we can utilise the spatial 
dependence of the energy level structure of the optical 
transitions to achieve qubit frequency selective reset and 
initialisation, tasks which are essential for the ultimate 
fault tolerant operation of a scalable device. 
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Mean Qubit Resonance Frequency (GHz) 


34.3 


Neighboring Qubit Resonance Frequency Separation (MHz) 


280 


Motional Sideband Extent (MHz) 


3 


Neighboring Qubit Optical Readout Frequency Separation (MHz) 


187 



TABLE I: Relevant parameters for ten Yb + ions, with Bo = 
IT, b = lOOOT/m, v x = 1MHz, and h = 10/im, held in 
individually tailored microtraps for uniform separation and 
J— couplings. 



(b) 



FIG. 1: Qubit coupling strengths, Jij /2n, in Hz, for a ten 
Yb + ion string in a magnetic field B(x) — Bo + bx, with 
Bo = IT, b = 1000 T/m, where (a) the ion string is held 
in a single harmonic trap of strength 1MHz, and (b) where 
each ion is held at a uniform spatial separation of 10/im, via 
individually tailored microtraps with an end trap strength of 
1MHz. 
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FIG. 2: Microtrap frequencies, Uj, as a function of ion-chain 
length, N, for a uniform ion spacing of 10/im, end-trap fre- 
quency of 1MHz and magnetic field gradient of lOOOT/m. Mi- 
crotrap frequencies are symmetrically distributed about the 
geometric center of the ion-chain, curve #k refers to the mi- 
crotraps j = k + l,N — k, in a chain of length N. 
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